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The unitarity condition unambigously requires the Orear region to appear in between the diffrac- 
tion cone at low transferred momenta and hard parton scattering regime at high transferred momenta 
in hadron elastic scattering. It originates from rescattering of the diffraction cone processes. It is 
shown that such region has been observed in the differential cross section of the elastic pp-scattering 
at ^/s=7 TeV. The Orear region is described by exponential decrease with the scattering angle and 
imposed on it damped oscillations. They explain the steepening at the end of the diffraction cone 
as well as the dip and the subsequent maximum observed in TOTEM data. The failure of several 
models to describe the data in this region can be understood as improper account of the unitarity 
condition. It is shown that the real part of the amplitude can be as large as the imaginary part in 
this region. The overlap function is calculated and shown to be small outside the diffraction peak. 
Its negative sign there indicates the important role of phases in the amplitudes of inelastic processes. 



I. INTRODUCTION 

The TOTEM collaboration has published Q exper- 
imental results on the differential cross section of the 
elastic pp-scattering at the total cms energy ,/s=7 TeV. 
Among the most interesting features they observe the 
steepening of the diffraction cone near the squared trans- 
ferred momentum 0.3 GeV 2 , the dip at 0.53 GeV 2 and 
the maximum at 0.7 GeV 2 . We explain them as result- 
ing from the rigorous requirements of the unitarity con- 
dition. It prescribes the Orear regime characterized by 
exponential decrease with the scattering angle to start 
at transferred momenta just above the diffraction cone. 
The damped oscillations imposed on it lead to the dip 
in the differential cross section. No particular model has 
been used. 

At the same time there exist several models mostly 
based on reggeon approach. Their predictions are exten- 
sively cited in [lj. Being rather successful in the diffrac- 
tion cone, they fail to describe the new data quantita- 
tively beyond the diffraction peak. This demonstrates 
that the unitarity condition is not properly accounted 
there in these models. Since then some other models 
were proposed 0, Q . 

At the end of 60s the very first experimental data on 
elastic pp- and irp- scattering were obtained at energies 
between 6.8 and 19.2 GeV in the laboratory system 0- 
6]. They showed that just after the diffraction cone which 
behaved as a Gaussian in the scattering angle there was 
observed the exponentially decreasing with the angle be- 
havior which was called as the Orear regime in the name 
of its investigator [5] . Some indications on the shoulder 
appearing at the beginning of this region (evolved later 
to the dip at the higher ISR energies) were also obtained. 
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The special session was devoted to these findings at the 
1968 Rochester conference in Wien. 

The theoretical indications on the possibility of such 
regime were obtained even earlier (7H9| but the results 
did not fit new experimental findings. 

At the same time the simple theoretical explanation 
based on rigorous model-independent consequences of 
the unitarity condition was proposed [13, HH an d a care- 
ful fit to experimental data showed good quantitative 
agreement with experiment [l2l |. 

We follow these ideas to demonstrate that they are also 
applicable to the recent data of the TOTEM collabora- 
tion at the LHC at energies as high as 7 TeV. 



II. THEORETICAL DESCRIPTION 

The elastic scattering proceeds mostly at small angles. 
The diffraction peak has a Gaussian shape in the scat- 
tering angles or exponentially decreasing as the function 
of the transferred momentum squared 




where the four-momentum transfer squared is 

t = -2p 2 (l-cos0)«-p 2 2 (0<0 d <l) (2) 

with p and denoting the momentum and the scattering 
angle in the center of mass system and B known as the 
diffraction slope. 

At large energies the forward scattering amplitude has 
a small real part as known from the dispersion relations 
[13l Il4j. Therefore, in the first approximation, it is rea- 
sonable to assume that its real part is negligible within 
the diffraction peak < 0^. Then the elastic scattering 
in this region can be described by the amplitude 

A(p,8)^4zp 2 * t e- B P 2e2 / 2 (3) 
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with a proper optical theorem normalization to the total 
cross section a t in the forward direction. We stress that 



Eq. (|3]) follows directly from experimental results and 
does not appeal to any particular model. 

Let us have a look at the unitarity condition which is 



ImA(p,6) = I 2 (p,0)+F(p,0) 
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+ F(p,9). (4) 



The region of integration in is given by the conditions 

\0i-(h\<0, 9 < e 1 + e 2 < 2tt - e. (5) 

The integral term represents the two-particle interme- 
diate states of the incoming particles. The function 
F(p, 9) represents the shadowing contribution of the in- 
elastic processes to the elastic scattering amplitude. Fol- 
lowing Van Hove 0] it is called the overlap function. It 
determines the shape of the diffraction peak and is com- 
pletely non-perturbative. Only some phenomenological 
models pretend to describe it (see also [l5[ where its 
shape is obtained using the unitarity relation in com- 
bination with experimental data). 

Now, let us consider the integral term I 2 outside the 
diffraction peak. Because of the sharp fall-off of the 
amplitude §3§ with angle, the principal contribution to 
the integral arises from a narrow region near the line 
9 1 + 9 2 ~ 9. Therefore one of the amplitudes should be 
inserted at small angles within the cone while another 
one is kept at angles outside it. At the beginning, let 
us neglect the real parts of the amplitude both in the 
diffraction region and at large angles. We insert Eq. fijfy 
for one of the amplitudes in I 2 and integrate over one of 
the angles. Then the linear integral equation is obtained: 



lmA(p, 9) = 
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+F(p,d). (6) 



It can be solved analytically (for more details see |10llll|) 
with the assumption that the role of the overlap func- 
tion F(p, 9) is negligible outside the diffraction coneflfjj. 
To account for the real part of the amplitude, one re- 
places at by Otfp where f p = 1 + pdPi with average 
values of ratios of real to imaginary parts of the am- 
plitude in and outside the diffraction cone denoted as 
Pd and pi correspondingly. It follows from Eq. Q that 
A X A* 2 -> ImAilmA 2 (l + p x p 2 ). 

Using the Fourier transformation one gets the solution 



lmA(p, 9) = C e V °%fp y 
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This shape has been obtained from contributions due 
to the pole on the real axis and a set of the pairs of 
complex conjugated poles. Correspondingly, it contains 
the exponentially decreasing with 6 (or \f\t\) term (Orear 
regime!) with imposed on it damped oscillations. Let us 
mention the papers [l7| where non-damped oscillations 
were predicted in the reggeon exchange model but they 
are not observed in experiment. 

The elastic scattering differential cross section outside 
the diffraction cone (in the Orear regime region) is 



da 
P\dt 
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The first (Orear) term is exponentially decreasing with 9 
(or \/jt[) and the second term demonstrates the damped 
(n = I) oscillations which are in charge of the dip- 
maximum structure near the diffraction cone. The omit- 
ted terms with larger n in Eq. are damped stronger 
because they contain ^/n in exponents. Let us note that 
the exponents of the damped terms are much larger nu- 
merically than that of the Orear term if the experimen- 
tally measured values of the diffraction cone slope B and 
the total cross section a t are inserted. Namely B and 
a t determine mostly the shape of the elastic differential 
cross section in the Orear region between the diffraction 
peak and the large angle parton scattering. The value of 
AirB/at is so close to I that the first term is very sensi- 
tive to pi . Thus it becomes possible for the first time to 
estimate the ratio pi from fits of experimental data. 

Beside the overall normalization constant p\ this for- 
mula contains the constants p 2 and <p which determine 
the strength and the phase of the oscillation [l8j]. They 
can be found from fits of experimental data. The con- 
stant p\ is determined by the transition point from the 
diffraction cone to the Orear regime. The constants p 2 
and 4> define the depth of the dip and its position. 

What concerns the ratios p's, one can choose pd ~ 0.14 
as prescribed by the dispersion relations for its value at 
t = [HI, [3| and use pi as another fitted parameter 
which influences the exponents in Eq. ([9]). 

Let us note that all parameters can depend on energy 
as well as the values of the diffraction cone slope B and 
the total cross section cr t . Surely, this is unimportant if 
the fit is done at a fixed energy as in the present paper. 

The unitarity condition is not a complete theory. It 
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imposes some restrictions on its consequences however. 
Its solution predicts the dependence on p9 ^J\t\ but 
not the dependence on the collision energy! Nevertheless, 
main exponents in ([§]) depend on energy. We are able to 
predict them at different energies if the dependence of 
the diffraction slope B and the total cross section a t is 
known from experiment. In this way different reactions 

I 



(including pp, in particular) may be analysed. 

Apart from comparison of theoretical predictions with 
experimental data one can get some knowledge about 
the overlap function F(p, 9) (see [15]). It is important, in 
particular, to confirm the assumption about its smallness 
outside the diffraction peak. Then the equation Q is 
used as an expression for F(p, 9): 



F(p, 9) = 16p 2 



1/2 



8/(1 + p d p) 
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da da 
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(10) 



where z, = cos ( 



K(z,z 1 ,z 2 ) = l — z 2 — z 2 —z 2 +2zz 1 Z2, 
and the integration limits are zf — zz-i ± [(1 — z 2 ){l — 

zl)\ 1/2 - 

At y/s = 7 TeV the angles are extremely small so that 
the kernel becomes very singular. K is close to but 
integrable. The divergence is of the type J dzjsfz and 
can be computed. Computing F in the diffraction cone 
one uses p — p c i- Outside it p = p\. 

Let us mention that the inhomogeneous equation (|4]) 
has been solved [ll[ by iterations with the overlap func- 
tion approximated by F/sai ne i = cxp[— Bi n p 2 9 2 /2]. Its 
more precise approximation is required to get accurate 
results but it is important that the conclusion about the 
phase <fi remains valid. 

Below we show and discuss the obtained results. 



III. A FIT OF THE EXPERIMENTAL DATA 
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FIG. 1. Fit of experimental distribution of elastic pp- 
scattering at y/s — 7 TeV 



Having at our disposal Eq. © we try to fit exper- 
imental distribution of elastic pp-scattering at y/s = 7 
TeV. The experimental values of B = 20.1 GeV -2 and 
a t = 98.3 mb were used in ©. We expect that Eq. © 
must be applicable fom the end of the diffraction cone at 
\t\ = 0.3 GeV 2 to the beginning of hard parton processes 
at |t| > 1 GeV 2 . The result is shown in Fig.[TJ 

It is seen that the fit is quite successful in the expected 
applicability region. First of all, we notice the steeper de- 
crease in the region 0.3 < |t| < 0.36 GeV 2 compared to 
the slope of the diffraction cone at \t\ < 0.3 GeV 2 as 
observed in experiment. It is explained here as the nega- 
tive contribution of the oscillating term in Eq. ([9]) . That 
determines the phase <f). The dip develops at \t\ — 0.53 
GeV 2 where the cosine in the second term is close to -1. 
Then this term increases, becomes positive and leads to 
the maximum at \t\ « 0.7 GeV 2 . The positions of the 
dip and of the subsequent maximum are uniquely deter- 
mined by the period of oscillations At = 2ir/B which is 
predicted by the unitarity condition and depends only on 
the well measured slope of the diffraction peak B. The 
damping exponent in front of the cos-term becomes so 
strong at larger |t| that the simple Orear regime with 
the first term in Eq. © prevails. Let us note that the 



exponent in this term is rather small because the ratio 
inB/at is very close to l[TJ[. Therefore it is extremely 
sensitive to the parameter pi. That helps determine this 
parameter. 

Hardly any oscillations will be observed at large \t\. 
The exponent in the oscillation term is very large and 
strongly damps it. One could pretend to observe the 
next weak oscillation at \t\ ~ 0.9 - 1.0 GeV 2 . However it 
would require very high precision. It is interesting to note 
that the damping increases with energy due to increase of 
the slope B. At the same time the shrinkage of the cone 
leads to the shift of the Orear regime (and the dip) to 
smaller angles at higher energies so that the oscillations 
are stll noticeable there. 

Let us list and discuss the parameters in Eq. (j9]) which 
we found by the fitting procedure: p± = 18.71; pi = 
115.6; = -0.845; pi w -2. The large value of 
pi demonstrates that the dip is well pronounced in the 
data. Up to now the only possible mo del- independent 
estimate of the ratio of real to imaginary parts of the 
elastic scattering amplitude was available from the dis- 
persion relations at t = 0. It is for the first time that it is 
done at large \t\ in a model-independent way and shows 
that this ratio is of the order of 1 there. Surely, there 
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are many models where this ratio is calculated in a wide 
range of i-values. There is no common consensus about 
their validity however. The parameter (j> is so close to its 
theoretical estimate that it was not even necessary to use 
it as a free one. 

Now we discuss the role of the parameters. 

1 . The parameter p\ is in charge of the overall normal- 
ization and, consequently, of the smooth transition 
from the diffraction cone to the Orear region. 

2. The parameter p2 defines the amplitude of the os- 
cillations and, consequently, the depth of the dip. 
In combination with <fr it leads to the steepened 
slope at 0.3 < \t\ < 0.36 GeV 2 . 

3. The phase <f> determines the position of the dip and 
the beginning of the steepened slope. Actually, it 
was shown in [ll| that it can be obtained from the 
iterative solution of the non- linear equation It 
is almost independent of the form of F(p, 9) so that 
\4>\ ~ 7r/4. Nevertheless this problem asks for 
further studies. 

4. The parameter pi in f„ is in charge of the exponen- 
tial slope at |t| above the maximum (together with 
B and at in the first term of ©)■ It is negative 
and rather large (in the absolute value). 

5. The relative position of the dip and the maximum 
(the period of oscillations) is determined only by 
the diffraction cone slope B (the second term in 
©)• 

IV. THE OVERLAP FUNCTION 

As follows from experiment, the inelastic cross section 
is much larger than the cross section of elastic scattering 
at high energies. Therefore the overlap function is much 
larger than the integral term in the unitarity relation 
at small t. To see what is the contribution of inelastic 
processes to the unitarity relation at any values of t, it is 
instructive to calculate the overlap function according to 
Eq. (fT0|) . In [15J that has been done at the assumption of 
ratios of real to imaginary parts p equal to zero both at 
small and large t. Now with the above estimate of pi we 
can take it into account. Nevertheless, the calculations 
were done with and without account of p to compare with 
previous results and estimate the role of p. The results 
are shown in Fig. [21 

There are several distinctive features observed. First 
of all, as expected, the overlap function drops down very 
fast with increase of the transferred momentum \t\ and 
determines the shape of the diffraction cone. Second, 
it crosses the abscissa axis at \t\ rj 0.3 and becomes 
negative. Namely there the Orear regime starts work- 
ing. If compared to low energies [l5j], the overlap func- 
tion becomes narrower at higher energies. Third, it is 
small and changes very slowly outside the diffraction cone 
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FIG. 2. The overlap functions calculated with pa — pi — 
and with pd = 0.14; pi — —2 (closest to the abscissa axis). 



similarly to the low energy behavior. Intuitively, this 
smallness may be understood as a consequence of strong 
destructive interference between amplitudes of inelastic 
processes with very different kinematics. In one of these 
amplitudes the final state must be turned to the large an- 
gle 6 relative to the direction of initial particles. Thus the 
overlap of these two processes is small. Fourth, the ac- 
count of p does not change qualitatively this conclusion in 
general even though somewhat changes the numerical es- 
timates diminishing |.F| further. This follows from better 
fit of experimental data with pi different from zero. Fifth, 
the negative sign of F imposes a severe problem to the- 
orists because it shows the important role of the phases 
of matrix elements of inelastic processes and their strong 
interference when trying to reconstruct elastic scattering 
from two inelastic processes turned by t one to another. 



V. CONCLUSIONS 

Thus we conclude: 

• At intermediate angles between the diffraction cone 
and hard parton scattering region the unitarity con- 
dition predicts the Orear regime with exponential 
decrease in angles and imposed on it damped os- 
cillations. Earlier, this solution was helpful in ex- 
plaining this regime at lab. energies 8-20 GeV. 

• The experimental data on elastic pp differential 
cross section at y/s=7 TeV in this region are fit- 
ted by it with well described position of the dip at 
\t\ w0.53 GeV 2 , its depth and subsequent damped 
oscillations with the predicted period about 0.3 
GeV 2 . The large amplitude of the oscillations and 
their negative sign explain the steepened slope at 
0.3 < \t\ <0.36 GeV 2 . The positive sign of the os- 
cillating term at \t\ ~ 0.7 GeV 2 leads to the maxi- 
mum. Strong damping of the oscillations at higher 
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values of \t\ results in clear signature of the sim- 
ple exponential (in y/\t\) behavior observed first by 
Orear which extends up to \t\ as 1.5 GeV 2 . 

• A good fit allows without using any definite model 
for the first time to estimate the ratio of real to 
imaginary parts of the elastic scattering amplitude 
in this region (pi rj —2) far from forward direction 
t=0. 

• The overlap function at 7 TeV has been calculated 
using only the experimental differential cross sec- 
tion and the above estimate of the ratio of real to 



imaginary parts. As at low energies, it is small and 
negative in the Orear region. That confirms the as- 
sumption used when solving the unitarity equation 
and shows that the phases of inelastic amplitudes 
become crucial in any model of inelastic processes. 
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